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This paper estimates an unobserved components model to explore coal consumption in the USA. We 
ask whether coal consumption exhibits hysteresis, defined as a dynamic structure in which the cyclical 
component of coal consumption has permanent effects on the natural component. In contrast to pre- 
vious analysis, we provide evidence in favor of hysteresis in coal consumption, by using the nonlinear 
framework proposed recently by Pérez-Alonso and Di Sanzo [13], in which threshold type nonlinearities 
are introduced by allowing past cyclical consumption to have a different impact on the natural compo- 
nent depending on the regime. The article discusses implications of the findings for energy conservation 
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1. Introduction 


Efforts to find alternative and cleaner energy sources are inten- 
sifying given the public interest in the reduction of the negative 
effects of carbon dioxide emissions. But meanwhile governments 
around the world regard energy conservation policies as a promis- 
ing candidate in this respect. In this context, the study of the effects 
of these conservation policies on the reduction of energy consump- 
tion is a highly hot topical issue at the time of writing, given that 
these interventions often impose sizeable costs on the taxpayer. 

Given these costs, the lack of robust evidence associating these 
policies with long term effects on the consumption is particularly 
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striking.! At the heart of this question is how energy consumption 
evolves. If energy consumption is trend-stationary, policy shocks 
such us the energy conservation policy can be regarded as tran- 
sitory: energy consumption eventually reverts to its underlying, 
long-run (“natural”) component. For instance, one can conceive 
of an energy conservation policy, which thereby creates a down- 
ward shift in consumption. If energy consumption is stationary, 
the shock dies away once the policy has been implemented, and 
energy consumption settles up at its new, lower, level. If on the 
other hand the energy consumption is non-stationary, such shocks 
can have permanent effects. For example, the energy conservation 


1 The available evidence is mixed. Some authors suggests that energy consump- 
tion is stationary in which case policies oriented toward the reduction of energy 
consumption only have short-term effects [4-7] while others [8,9] provides evi- 
dence on the contrary. 
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policy considered above, decreases energy consumption in a way 
which affects all subsequent levels of consumption. A necessary 
(but not sufficient) condition for this to occur is a unit root in energy 
consumption. 

In a time-series context, hysteresis can be defined and measured 
in various ways. A popular approach in the empirical literature 
simply equates hysteresis with the existence of a unit root in a 
variable (see, Rged [1], for a survey). An alternative approach pro- 
posed by Jaeger and Parkinson [2,3] posits a more demanding 
criterion: hysteresis exists if shocks (such as the one-off policy 
change discussed above) affect the natural component of a variable, 
which itself follows a unit root process. In which case, temporary 
shocks have permanent effects while the cycle does not evolve 
independently of the natural component; it then follows that a 
unit root is a necessary but not a sufficient condition for hystere- 
sis. In this article, we adopt Jaeger and Parkinson’s [2,3] definition 
of hysteresis in order to conduct a searching test and to explore 
whether energy conservation policies on coal consumption have 
long-term effects. 

To test for hysteresis in this way, we decompose coal con- 
sumption into two unobservable components: a non-stationary 
“natural or permanent” component, and a stationary “cyclical” 
component. These components can be estimated by maximum like- 
lihood using the Kalman filter. Although Jaeger and Parkinson’s 
approach has been applied extensively in labor economics [10-13] 
to the best of our knowledge its application to energy economics 
is novel. 

Once hysteresis is tested in the linear model, the new test for 
hysteresis based on a nonlinear unobserved components model, 
proposed by Pérez-Alonso and Di Sanzo [13], is applied. This test 
introduce, threshold type nonlinearities by allowing past cyclical 
component to have a different impact on the natural component 
depending on the regime of the economy. Note that we must take 
into account this possibility because the estimation of linear rela- 
tions could yield spurious inference results if the nonlinear model 
provides a better empirical description of our data. 

A further motivation for our empirical analysis is that hysteresis 
exists this implies that energy conservation policies might be more 
powerful than has been thought hitherto. 

This article has the following structure. The next section 
describes the data and the estimation methodology. The third 
section presents and discusses the results. The final section con- 
cludes with a discussion of policy implications and some promising 
avenues for future research. 


2. Data and methodology 
2.1. Data 


The data used are quarterly observations from 1973:1 to 2010:3. 
The coal consumption (measured in thousand short tons) and GDP 
data (measured in billions of chained 2005 dollars) are extracted 
from the U.S. Energy Information Administration (EIA) and the U.S. 
Bureau of Economic Analysis (BEA), respectively. Before conducting 
the empirical analysis data were seasonally adjusted. 


2.2. Econometric methodology 


As we mentioned before, several empirical studies equate hys- 
teresis in a time series with a unit root process. Others argue that 
hysteresis arises when changes to the cyclical component of a time 
series, ae induce permanent changes in the “natural” component 
of the series, CN. This is different to a unit root process. To com- 
prehend the different estimation strategies these approaches call 
for, decompose the series C; into the sum of its two (unobservable) 


components: the non-stationary natural component, CN, and the 
stationary cyclical component, ce : 


Ce = CN + CE (1) 


Now define the natural component as arandom walk plus aterm 
capturing a possible hysteresis effect: 


CN =CN, + ACL, +e (2) 


where the £ coefficient measures, in percentage points, how much 
the natural component increases if the consumption experiences a 
cyclical increase of 1%. Evidently a unit root in the consumption C; is 
necessary but not sufficient condition for the existence of hysteresis 
since a unit root could be generated by an accumulation of shocks 
to the natural component C™ while at the same time £=0 [1]. In 
contrast, there is hysteresis if 6 >0. 

The specification of the model is completed by writing the 
cyclical component as a stationary second-order autoregressive 
process?: 


CE = pCO, + CO, + ef (3) 


where ¢ and p2 provide a measure of the periodicity of the cyclical 
component. 

To identify the model, the system is completed by augmenting it 
with an equation, which relates the cyclical component of the coal 
consumption and output growth, 


Dt = aD,_1 + CE + eP (4) 


where D; stands for the output growth rate at date t.3 

The random shocks eN, ef and eP are assumed to be mean- 
zero draws from the normal distribution with variance—covariance 
matrix 2; the state-space form of the model can be written as 


+(2)aa+($) (5) 
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To summarize, hysteresis is inferred if the coefficient £ is sig- 
nificantly different from zero. The coefficients of the model (4)-(7) 
are estimated by maximum likelihood using a Kalman filter. 

A non-linear version of this model (5)-(7) can also be estimated, 
to take account of the possibility that coal consumption displays 
asymmetries in adjustment dynamics in response to positive and 
negative shocks. Relaxing the linearity assumption may allow a 
better estimation of hysteresis if it exists. When we talk about “pos- 
itive” or “negative” shocks, we do so relative to some threshold 
level, t (where T is not necessarily zero). To explore whether asym- 
metries exist, we estimate a non-linear version of the unobserved 
components model by allowing past cyclical consumption to have 
a different impact on the natural consumption, which depends on 


2 The assumption of a purely autoregressive process for the cyclical equation can 
be relaxed in favor of more general (and possibly more parsimonious) autoregressive 
moving-average specifications. In the present application, an AR(2) fits the data best 
according to AIC comparisons. Full results are available from the authors on request. 

3 However, as Proietti [14] states, this model is identified without this additional 
equation. 
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the regime of the economy. Specifically, we replace the state-space 
equation (6) with the Threshold Auto Regressive (TAR) specification 


Py /1 0 0) (Ba) (BY B 
C |=|0 G1 p cC J+] 0 Ecl 0 Jie 
S 01 0 Cs 0 0 


where If and I; are the Heaviside indicator functions such that: 


or of 


This model can be estimated via maximum likelihood using the 
Kalman filter, where t is unknown so that it is estimated along with 
the other parameters of the model 6* and f-. In this context a test 
for asymmetry becomes a test for linearity, i.e. a test for a single 
regime against the alternative of two regimes. The null hypothesis 
we are interested in is Ho : B* = 6- vs. Ho: B* + B-.If we reject the 
null of linearity, there is evidence for the presence of a type of non- 
linear hysteresis in coal consumption, i.e. with cyclical shocks being 
propagated asymmetrically to the natural component. Given our 
model, the asymptotic distribution of conventional test statistics is 
not x2.4 To circumvent this problem we follow [13] who suggest 
using bootstrap methods to approximate the sampling distribution 
of the test statistic. 


3. Results 


This section presents the results in several stages. First, we test 
what Jaeger and Parkinson [2,3] have characterized as a necessary 
but not sufficient condition for hysteresis, namely the existence 
of a unit root in the coal consumption time-series. Second, we 
estimate the linear unobserved components model outlined in Sec- 
tion 2.2, incorporating a unit root as a maintained hypothesis. The 
third subsection explores the possibility of asymmetric behavior 
in adjustment dynamics, by estimating the nonlinear TAR unob- 
served components model, using the estimation strategy proposed 
by Pérez-Alonso and Di Sanzo [13]. 


3.1. Unit root tests 


In order to test the hypothesis of non-stationarity, we apply 
the traditional Augmented Dickey—Fuller (ADF) test and a modified 
version of the Dickey—Fuller and Phillips—Perron tests proposed by 
Ng-Perron [17]. This comprises a class of modified tests, M, with 
GLS de-trending of the data and use of the modified Akaike infor- 
mation Criteria to select the autoregressive truncation lag. Table 1 
reports the results of Ng—Perron tests, MZ¢4 and MZE“, originally 
developed in Stock [18] with GLS de-trending of the data as pro- 
posed by Elliot et al. [19]. In addition, Ng—Perron proposed a similar 
procedure that corrects the problem associated with the standard 
Augmented Dickey Fuller test, MSB@5 and MPT“, All test statistics 
formally examine the unit root null hypothesis against the alterna- 
tive of stationarity. 

The results in Table 1 show that the null hypothesis of non- 
stationarity cannot be rejected, regardless of the test. However, it 
is well known that structural breaks in time-series can lead to spuri- 
ous inferences ofa unit root. To deal with this possibility, we employ 
the Zivot and Andrews [20] minimum ADF-t(min-t) procedure. The 
min-t statistics reported in Table 2 show that the null hypothe- 
sis of a unit root in the time series still cannot be rejected. This 


4 See, e.g. Hansen [15] and Lo and Zivot [16]. 


buttresses our conclusion that a unit root exists in coal consump- 
tion. As noted above, a unit root is a maintained assumption 


N 
Et 


+ù @ (8) 


needed to test for Jaeger and Parkinson’s notion of hysteresis. We 
investigate the notion of hysteresis next. 


Table 1 
Unit root tests. 


Variable Coal consumption 

MZSS 0.492 

MZ% 0.519 

MSB 1.056 

MPT 69.198 

Lag length 11 

ADF —2.588"" 

Lag length 11 

Range 1973:1-2010:3 

Critical values (%) MZSS MZSS MSBeS MPT 
1 —13.800 —2.580 0.174 1.780 
5 —8.100 —1.980 0.233 3.170 

10 —5.700 —1.620 0.275 4.450 

Critical values (%) ADF 
1 —3.476 
5 —2.881 

10 —2.577 


Notes: Test statistics defined in the text. “Lag length” refers to the lag length used 
in the Ng-Perron and ADF tests, respectively. The critical values are tabulated in Ng 
and Perron [17]. 
*Rejects null hypothesis at 10% significance level. 
*“Rejects null hypothesis at 5% significance level. 

*“ Rejects null hypothesis at 1% significance level. 


Table 2 
Unit root tests allowing for structural breaks. 


Variable Coal consumption 

A Min-t —0.807 (1983:3) 
Lag length 12 

B Min-t —1.642 (2005:1) 
Lag length 12 

C Min-t —1.571 (2000:2) 
Lag length 12 

Range 1973:1-2010:3 


Notes: Periods corresponding to min-t statistics are indicated in parentheses. Critical 
values for the min-t are given by Zivot and Andrews [20]. Asterisks are as in Table 1. 
Min t-statistics are computed using sequential regressions over 1 < trend break<T 
based on the following equations: 


k 
Ax; = 644 ht + 8DU + ox + Soo Ax, j ter (A) 
j=1 
k 
Ax; = 68 + 68t + BDT + aPxii + Scop ax; +e (B) 
jel 
k 
Ax; = 8E +8°t + SEDU + 8SDT +.0°xt1 + S “af Ax; +e (©) 


jel 


where the dummy variables DU; =1 and DT, =t — TB for t> TB and 0 otherwise, and 
TB denotes the period at which a possible trend break occurs. Critical values for the 
min-t are given by Zivot and Andrews [20]. In model (A): 1% (—5.34) 5% (—4.80) 10% 
(—4.58); model (B): 1% (—4.93) 5% (—4.42) 10% (—4.11); model (C): 1% (—5.57) 5% 
(—5.08) 10% (—4.82). 
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Table 3 
Estimates of the linear unobserved component model. 


Natural rate equation 


B 0.000 (0.014) 

ON 0.021" (0.003 
Cyclical rate equation 

Qı 0.488" (0.291) 

Q2 —0.059 (0.071) 

oc 0.015%" (0.004 
Identification equation 

a 0.654" (0.090 

ô —12.466' (6.471) 

op 0.773" (0.092 
Range 1973:1-2010:3 


Notes: Standard errors are in parentheses. 
* Rejects null hypothesis at 10% significance level. 
**Rejects null hypothesis at 5% significance level. 
* Rejects null hypothesis at 1% significance level. 


Table 4 
Non-linear model estimation results. 


Natural rate equation 


B 1.076" (0.123) 2.028" (0.679) 
ON 0.029"** (0.003) 
Cyclical rate equation 
Qı 0.269” (0.130) 
Q2 0.407*" (0.104) 
oc 0.010” (0.004) 
Identification equation 
a 0.300" (0.159) 
ô 7.185™ (2.206) 
op 0.479” (0.211) 
Threshold 0.0137 
Delay lag 3 
% obs. 56.8 73.2 
Range 1973:1-2010:3 


Notes: Standard errors are in parentheses. 
* Rejects null hypothesis at 10% significance level. 
* Rejects null hypothesis at 5% significance level. 
* Rejects null hypothesis at 1% significance level. 


3.2. The linear unobserved component model 


Table 3 presents the results of estimating (5)-(7) for coal con- 
sumption. The parameter £ is zero and not statistically significant. 
However, should be rejected the existence of hysteresis in coal con- 
sumption? A likely reason of this rejection could be the presence 
of nonlinearity, and we should check it. 


3.3. Asymmetries 


We next check whether our results are robust to the linear speci- 
fication of the unobserved component model.’ This involves jointly 
estimating the structure (4), (6) and (7) to determine whether there 
is a threshold for cyclical consumption which is associated with 
asymmetric hysteresis responses. We wish to check whether the 
findings in the previous subsection are robust to possible asym- 
metries, or whether they were merely an artifact of the restrictive 
technical assumption of linearity. 

The p-values calculated following the bootstrap technique 
described in detail in Pérez-Alonso and Di Sanzo [13]. The null 
hypothesis Hp: 6*=8- is rejected. Then, the hysteresis effect is 
significant at the 1% level. 


5 In this case, the assumption of stationarity should be tested by using an alterna- 
tive method. We employed the Caner and Hansen [21] methodology to test for a unit 
root ina TAR model. The null hypothesis is that there is not a unit root. The p-value - 
calculated by using the bootstrap technique described in Pérez-Alonso and Di Sanzo 
[13] - is 0.04. 


In Table 4, we report the ML estimates of the hysteresis param- 
eters. The threshold parameter is t=0.0137. Hence, the threshold 
model splits the regression into two regimes depending on whether 
or not the threshold variable is higher than this threshold param- 
eter. In both regimes, the hysteresis parameter is statistically 
significant and positive and the one associated with Regime 2 
is greater than that of Regime 1. This result points to asym- 
metric responses of the natural component as regards cyclical 
consumption. Therefore, energy conservation policies might be 
more powerful than has been thought before. That is because any 
decrease in coal consumption brought about these policies will be 
incorporated into all future levels of consumption. 


4. Conclusions 


This paper estimated unobserved components models for coal 
consumption in the United States. Defining hysteresis in terms of 
the interdependent evolution of a non-stationary natural compo- 
nent and a stationary cyclical component, thereby distinguishing 
hysteresis from natural component shocks, the results provide 
robust evidence of hysteresis, although the hysteresis’ coefficient 
is time varying. This implies that policy shocks in have perma- 
nent effects on coal consumption, although with different intensity 
depending on the regime. Future work might fruitfully apply the 
methodology used in this article to a broader range of energy 
sources and revisiting the robustness of some previous studies on 
energy consumption. 
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